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Об устойчивости в моделях математической 
зпидемиологии на основе функций Ляпунова-Вольтера 


Робота присвячена побудові функцій Ляпунова-Вольтера в задачах математичної епідеміології. 
Розглядаються моделі 5ІВ, 51АВ та модель співіснування двох штамів вірусу. 
Ключові слова:математична епідеміологія, 51А-модель, функція Ляпунова-Вольтера. 


ТЬе ууогік 15 деуоїсд 10 пе соп5ігисіїоп ої Гуарипоу-Уойега Тапсіїопя іп ргобіешіє ої плаїпетанса! 
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Работа посвящена построєнию функций Ляпунова-Вольтера в задачах математической зпидемиологий. 
Рассматриваются модели 51КВ, 51АК и модель сосуществования двух штаммов вируса. 
Ключевье слова: математическая зпидемиология, 5ІВ-модель, функция Ляпунова-Вольтера. 


Другий метод Ляпунова є одним із найконструктивніших підходів до дослідження 
стійкості динамічних систем. Перевагою методу є те, що в певних випадках систем 
(наприклад, лінійних) вдається отримати не лише достатні, але й необхідні умови 
СТІЙКОСТІ. 

У той же час проблемою залишається конструктивна побудова функцій Ляпу- 
нова та отримання умов стійкості для конкретних динамічних систем. 

Одним із важливих застосувань є задачі математичної епідеміології на основі 
518 -моделювання. Слід зазначити, що проблеми стійкості таких систем вивчаються 
в роботі (1). При цьому традиційно використовується підхід дослідження стійкості за 
першим наближенням. Отримані умови стійкості при цьому виражаються в термінах 
показника відтворюваності. 
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Метою даної роботи є розробка підходу отримання достатніх умов стійкості на 
основі функції Ляпунова-Вольтера. Перевагою такого підходу є врахування нелі- 
нійного характеру моделі. 

Найдавнішу історію в моделях популяційної динаміки має функція Ляпунова 
«типу Вольтера»: 


у «У ах -хіпх,) 
ізі 


Дана функція застосовується в моделях Лотки-Вольтера |2-5| і вперше була 
запропонована самим Вольтером, хоч і він не вживав термін «функція Ляпунова». 

Функція Ляпунова-Вольтера почала активно впроваджуватися в епідеміо- 
логічних моделях, починаючи з 1980-х років. 

Завдяки роботі (6) функції Ляпунова-Вольтера почали використовуватися для 
вивчення стійкості у багатовимірних системах, що грунтуються на законі дії мас. 
Складністю залишається вибір коефіцієнтів, що гарантують від'ємновираженість 
похідної |7|. В даній роботі розвивається загальний метод вибору коефіцієнтів 
функції Ляпунова-Вольтера, що грунтується на нерівності між середнім арифме- 
тичним 1 геометричним. 

Функція Ляпунова в 5ІК-моделі. 

Розглянемо 5ЇВ-модель, представлену діаграмою рис. І. 


их В15 а 
и5 Ш ик 


Рисунок 1 - Функція Ляпунова в 5ІК-моделі 


Модель задається такою системою звичайних диференціальних рівнянь: 
5" (М - 5)- РІЗ, 
Г - ВІ5 -(шча)|, а) 
В'-аї - В. 
Із загальної теорії відомо, що існує єдиний ендемічний стан рівноваги, що за- 
довольняє ендемічні співвідношення: 


м зе ВГ5, 
ВГ8" «(шча)і?, 0) 
аГ зб". 


Розглянемо функцію, що претендує бути функцією Ляпунова: 
Ук - (5-51 5)-ЩІ - Г'ї 1) 

де 420 - поки що невідома стала. 

Її похідна на траєкторіях (1) має вигляд: 


пози 5)- ВІЗ шоб ния 


з-4|ВІЗ - (и за) - ВГ85-(шча)г| 
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Використовуючи ендемічні співвідношення (2), маємо: 
а щі ж ВГ8" и Й «(дя ве) из п б 


ж 


з-(а -1)рі5 ««Ї-(шчаді - ВІГ8 - ВГ'5"|- ша - Х| у 
ч яз за) 5 а 2 | з-(а -1)рІ8 -- ЦВУ -а(ш-ка)) 
З другої рівності в (2) бачимо, що 
85" «(ш ча) 


Отже, поклавши 4-1 і використавши нерівність між середніми арифметичними 
і геометричними бачимо, що Й, є від'ємною, тобто Й, 2 0. 


Розглянемо ЗІАВ-модель, представлену діаграмою рис. 2. 


ША 


Рисунок 2 - 51АВ-модель 


На її основі розглядається система диференціальних рівнянь: 
5" а М - 5)- ВИ - 64)5, 
І - ВІЗ -(ш- а), 


3 
4 - ВА5 -(и чт), У, 
В'-аї-пА-ив. 
Ендемічний стан рівноваги визначається ендемічними співвідношеннями: 
им з б ВГ 845, 
ВГ5 «(шча)ї, 


ВбА8" «(рчті), 
а зпА зшВ. 
Розглянемо функцію-претендент та функцію Ляпунова: 
Ух - (5-5 п5)-(1- Гл 1)-ЩА- А п А) 
де 420 - невідома система. 
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Її похідна та траєкторія (3) має вигляд: 
лан хо 5)- В(Г- 64)5 - нм в'я В(Т мі у 


- |рі5 -(шча) - ВГ5 «(и за) рзА8 -(нтп)А- ВбА"5 -(ш зл)4| 
Використовуючи ендемічні співвідношення (4) маємо: 


ух, я |рв" зо ВГ зн виг) з ші оо В би) - 


5 
-(р8' ВГ «ву Ро и" В(І бл) 

--|ВІ5 -(шчаді - ВГЗ" 5 вн 

м м -(н нп)л- рім'8-5 нрабз| - 2 2 "б 


ж ж 


кож 50 5 оч 5 5 
- ВІ/5 ан З уви 5 У є 15 в550 -4)А- 
- А(рі5" (шт) 
З другої рівності в (4) бачимо, що 85 «шча,аз третьої рівності в (4) 
Вб5 «шчтп. 
Отже, поклавши 4 - І і використавши нерівність між середніми арифметичними 1 
геометричними бачимо, що Й, є від'ємно визначеною, тобто Й, 2 0. 
Функція Ляпунова в моделі співіснування двох штамів вірусу. 
5 М - 5)- (ВИ, -81,)5, 
І «В51, «(р над)ї,, 
Кк заз (ш ї 0в.В, 
Уї- 0,В,К, І, З (ш Ея а, 


туг і, / -1,2,і є /. 
Із загальної теорії ми знаємо, що існує єдиний ендемічний стан рівноваги, який 
задовольняє рівняння: 


им з б" (ВІ ВГ, 
перо «(шана )/), 

о «(ша 0,8,Г)к, 
9.В.КГ -(рча,)ї,, 


і"3)7і 


(5) 


туг і, / -1,2,і є /. 
Розглянемо функцію Ляпунова: 
(1 -Глаг)ч 


І 


У» - (5-5"1п5)з) 


2 
ізі 


2 2 
зі 
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Її похідна вздовж траєкторій має вигляд: 


Й їв б мм -5)-(ВИ, орлом виб ВІ ря 
«УДряи- (ш на), - В75 (ша); їн 


КІ. 
З У мі чав - чі - Чинорадж 
Ї, 


і-1,іж) 


2 5 
«Умдлау - (шана, )ї, - у --8,8,8.1, -"(шча,) 
ізі 
Використовуючи ендемічні співвідношення (5) маємо: 


Ре Їр5 (ВГ кВ. 7)-ш5- (ВІ чн В5)8- 


Бе й (ша з (в зввк)є 5 -вВІ5 ува 
3 851, -(шча, ВК нина р 


- 

Ким 

8 

У 

з 
реа З 


К, гом я 5 
нач, В: дек) гг аноду жі 


ізн і і Про ооуе Ек б: В, Г 
- КІ, й є г; 
40 о указав нуво Пав У 
2, б В. Я В; Еби (ш-н а) У у В Іі ЕК і -(шча,) 
де (шана, аву, 
ЗЛ 
| І 
чн -(рч0,вГ)к, 
(чад уз ОВК, 
(ш та; у а 98 


ізі 


2 й 5 5" и тая 5" 5 
У сао нано ВВОДУ ЗЕ 


- Й та ува 
- 08 ГЕ С чо я ИН, с но В Я 
2 ,В; і р Г г / Ек І / у КІ 1 
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Враховуючи третє 1 четверте ендемічні співвідношення (5) маємо: 


. ж 2 ж 
2ц5112- 5 5 р Г8" р З 5 ка 
я 5" 5 -а 5 


п 57 
3 со 7 МН На У БК І, ВІ у, 
бо Брав р ЛЯ с ее ре 
жо ВО УВО ОБОКОБО ЛЬВОВ 


ГУ жо ц) 


Отже, в роботі вивчаються питання застосування другого методу Ляпунова до 


вивчення стійкості епідеміологічних моделей. Апробовано метод обгрунтування 
від'ємної визначеності похідної функції Ляпунова-Вольтера з використанням 
нерівності між середнім арифметичним і середнім геометричним. 
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ВЕ5СМЕ 
А.М. КисПпуага 


Оп Уабійу іп Моаєіз о! Матетаїіга! Ерідетіоіоєу Вазеа 


оп Гуарипо»- ИоПекга Кипсіїоп 

І єїмеп агбісіє деаї5 ут (Бе 5есопа тешодй ої Гуарипоу. ПП 15 опе ої Ше тобі 
сопзігисйуе арргоаспез Ко 5киду їБе з:аьШсу ої дупатиса! 5узіетя. ТПе адуапіаєе ої їБі58 
тефодй 15 "Ваї їп 50те са5е8 ої зуз5іеті5 ме сап єеї пос опіу 5иїНсіепі Би! а5о песез5агу 
сопаїнопя ої з:2ьШу. Тре іуевідануе арргоасії ої 5:абіШту ої їБе бт5і арргохітабоп 15 
ц5ед. Тре5е 5габійсу сопаїноп5 из ехргез5ей їп егтл5 ої Ре гаке ої гергодисібійсу. 

Опе шпрогіапі арріїсасіоп 15 "Бе ргобіет ої паїретагіса! ерідетіоїобу Ба5ей оп 
518.- тодейпє. Ії 5роцід Бе покей ШВаї Бе ргобіеті ої 5табїіу ої 5исП зузіетя аге 5їифтед 
їп біз гезеагсії Бу а соп5ігисіїує птлекпод ої 8киду ої пегайуєе деПпігепез8 ої Бе Фегіуанує 
ої Фе Гуарипоу-Уоїкетта їапебоп, ц5іпе бе іпедцайту Бебуееп Фе агіїптекіс пап апі 
Бе деогаєїтіс плеап. ТРе рарег 15 деуоїед іо їфе сопзігисійоп ої ПГуарипоу-УМоїегта 
Гипснопз5 їп ргобіетя ої плаїретакіса! ерідетіоіову. Тре тодеія 5ЇВ, 51АВ апа тодеї ої 
соехі5іепсе ої бо 5ітаїп5 ої Фе уїгиз аге ехатіпей. ТРре айт ої (Пі8 гезеагсі 15 їо демеїор 
ап арргоасі ої обіаїпіпе 5аййсіепі сопаїйопо ог 5:абійу Ба5ед оп Гуарипоу-УМоПетта 
Гипсноп. ТРре адуапіаєе ої'їБі8 арргоасі 15 Бе іпсогрогайоп ої поп-Ппеаг тоаеівз. 

У/е зкиду Ше дцезіоп ої арріїсайоп ої Ше 5есопа Гуарипоу тейод (о 5іиду їРе 
зіабійу ої ерідетіоіодіса! плодеіз. Аргобиуеїзуа ТРре теїрой ої 5їиду ої пебайує 
дейпіїепе55 ої Де Фегіуайує ої Ше Гуарипоу-Моїйега Гапсіоп ц5іпє їфе іпедцайу 
Бебмееп Фе агійптебіс плеап апа Бе беотеїгіс теап 15 Ке5іей. 


Стаття надійшла до редакції 23.12.2013. 
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